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Abstract 

We detected 2.8 bright points (BPs) per Mm^ in the Quiet Sun 
(QS) with the New Solar Telescope (NST) at Big Bear Solar Observa- 
tory; using the TiO 705.68 nm spectral line, at an angular resolution 
0.1" to obtain 30 min data sequence. Some BPs formed knots that 
were stable in time and influenced the properties of the granulation 
pattern around them. The observed granulation pattern within 3" 
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of knots presents smaller granules than those observed in a normal 
granulation pattern; i.e., around the knots a suppressed convection is 
detected. Observed BPs covered 5% of the solar surface and were not 
homogeneously distributed. BPs had an average size of 0.22", they 
were detectable for 4.28 min in average, and had an averaged contrast 
of 0.1% in the deep red TiO spectral line. 
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1 Introduction 

Observations of the solar photosphere 
performed in the G-band have revealed 
a plethora of tiny bright features, usu- 
ally concentrated in active regions or 
bordering the supergranules in the QS 
(Sanchez Almeida et al. 2001, Berger 
et al. 2004, Rimmele 2004, Sanchez 
Almeida et al. 2004, Beck et al. 2007, 
Utz et al. 2009, Viticchie et al. 2009, 
Sanchez Almeida et al. 2010, Viticchie 
et al. 2010). Near disk center, they ap- 
pear as "Bright Points" (BPs), roundish 
or elongated bright features located in 
the intergranular dark lanes (Dunn and 
Zirker 1973, Mehltretter 1974, Spruit 
1979, Title et al. 1987, Sanchez Almeida 
et al. 2004). G-Band BPs are known 
to be associated with strong magnetic 
flux concentrations of about 1.5 x 10^ 
G, hence, they can be used as a proxy to 
track small scale kG fields (Berger et al. 
1995, 1998, van Ballegooijen et al. 1998, 
De Pontieu 2002, Beck et al. 2007,Vitic- 
chie et al. 2009, Viticchie et al. 2010). 
Thus, investigation of the properties of 
BPs can shed new light on the abun- 
dance of kG fields in the quiet Sun. 

Simulations of magnetohydro dy- 
namics successfully reproduced the ob- 



servable effects of magneto-convection 
on the solar surface. Models predicted 
formation of small-scale concentrations 
of the magnetic flux caused by the inter- 
action of magnetic fields and convective 
motion (Spruit et al. 1990, Solanki, 
1993, Shelyag et al., 2004, 2007). MU- 
RaM simulation code (Vogler et al. 
2005) predict appearance of an increase 
in temperature in regions of down-flow 
(Shelyag et al. 2004 - Fig 3.) in the QS. 
A 3D LTE model (Hayek et al. 2010) 
provides a valid representation of the so- 
lar photosphere when scattering is not 
included in modeling of radiative trans- 
fer. Danilovic et al.(2010) demonstrated 
quantitative agreement of synthetic ob- 
servations made with MURaM code and 
real observations, indicating the valid- 
ity of the simulations and the inter- 
pretations of the spectro-polarimetric 
observations. 

Additionally, Morinaga et al. (2008) 
observed that magnetic flux concentra- 
tion may suppress convection in the ac- 
tive region; the strength of the magnetic 
flux is less relevant to the suppression 
of convection than the concentration of 
magnetic flux tubes. The BPs are the 
smallest flux tubes that can be investi- 
gated because their sizes are often com- 
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parable to the resolution of the instru- 
ment used. The size of the BPs varies 
from 0.17" (Berger et al. 2004) to 0.8" 
(Berger et al. 1995), while Beck et al. 
(2007) found that the effective diameter 
of BPs is 0.2", concluding that very few 
of the BPs are larger than 0.4". Utz et 
al. (2009) acquired two different results 
for the average diameter of BPs; 218 
km ( 0.3") when diffraction element was 
sampled at 0.108 arcsec/pixel and 166 
km ( 0.22") when sampled at 0.054 arc- 
sec/pixel. The authors emphasize the 
influence of resolution on the apparent 
size of BPs. Mehltretter (1974) found 
that BPs have typical dimensions of 100 
to 200 km (from 0.13" to 0.27"). 

Several methods for automatic de- 
tection of BPs have been developed 
(Bowelet and Wiehr 2001, 2007; Crock- 
ett et al. 2009; Utz et al. 2009). 
Sanchez Almeida et al. (2004) used 
a manual method to measure 0.3 BPs 
per Mm^ in the interior of the super- 
granulation cell, implying that BPs are 
ubiquitous in the photosphere. Sanchez 
Almeida et al. (2010) measured 0.97 
BPs per Mm^ using the same method 
but better resolution. The number BPs 
detected is highly dependent on the res- 
olution of the instrument and method 
used, consequently newer instruments 
and methods yield more BPs. 

In this paper, we present evidence of 

a mutual interaction between small scale 
bright features and granulation. We re- 



port the observational properties of BPs 
(i.e., size, velocity and contrast) and 
compare them with previous studies. 

2 Data and analysis 
methods 

Observations were performed with the 
New Solar Telescope (NST) at Big Bear 
Solar Observatory (Goode et al. 2010) 
on 29 July 2009. Photometry data of 
QS at the disk center were obtained us- 
ing an optical setup containing a TiO 
broadband filter and PCO.2000 camera 
(Cao et al. 2010) in the Nasmyth focus. 

Data consists of a sequence of 120 
bursts of 100 images. Each frame has 
an 10 ms exposure time with a tempo- 
ral cadence of 15 s between bursts. Im- 
ages have a sampling of 0.037" per pixel, 
which oversampled the diffraction ele- 
ment in the TiO spectral line. The field 
of view (FOV) encompassed 70" x 70" 
at the center of the solar disk. Image 
reconstruction and alignment of frames 
in the time series reduced the FOV to 
54" X 54". Despite excellent seeing con- 
ditions, we were not able to reach the 
diffraction limit of 0.1" for TiO 705.68 
nm spectral line, because NST is still in 
commissioning phase. 

Images were speckle reconstructed 
based on the speckle masking method 
(von der Liihe 1993). For this pur- 
pose we used the Kiepenheuer- Institut 



^Thc' software package offers a GUI interface to choose reconstniction method and set 
parameters of reconstruction. We used the Triple Correlation Method witli parameter 
settings: subfield = 5, Phase rec. limit=100%. Max. Iteration=30, SNR threshold=90, 
Weighing exponent=1.3 to complete image reconstruction. Wiener Noise Filter and Adap- 
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Speckle Interferometry Package (KISIP; 
Woger et al. 20080. 

After speckle reconstruction, images 
were aligned using a Fourier routine. 
This routine uses cross-correlation tech- 
niques and squared mean deviations to 
provide sub-pixel alignment accuracy. 
However, we did not implement sub- 
pixel image shifting to avoid substantial 
interpolation errors that sometimes ac- 
company this technique. 

We focused on the properties of 
small scale bright structures, e.g., their 
contrast, dimension, and lifetime. 

The term contrast is used as a mea- 
sure of relative brightness as compared 
to the intensity of the QS, and is defined 



as: 

C=^^, (1) 

where Im is the mean intensity of all pix- 
els from the targeted image area and Iq 
the mean intensity of the truly QS in 
our FOV. FOV contained a plethora of 
BPs so we defined QS in following way: 
As a QS we chose subareas that did not 
contain any of the resolvable small struc- 
tures (Fig{T]A). Intensity from these ar- 
eas was used to obtain the QS intensity. 
This provided us with the intensity of 
QS that had no resolvable contribution 
of the small structures. Fig. [TjB shows 
an example of one of the studied objects. 



BP knot 





Figure 1: Panel A presents an example of the quiet sun as we defined it for 
this work. Panel B presents one of the stable locations of the BPs in our field 
of view, as the contrast to the QS. 



tive low pass filters were used. 
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We measured diameter of BPs when 
they were the brightest to achieve 
greater accuracy and took intensity pro- 
files using full-width at half maximum 
(FWHM). For aspheric BPs we mea- 
sured FWHM from light profile on the 
longest axis. 

We calculated feature velocities us- 
ing the nonlinear affine velocity estima- 
tor (NAVE) method (Chae and Saku- 
rai, 2008), which applies an affine veloc- 
ity profile instead of a uniform velocity 
profile commonly used in the local cor- 
relation tracking (LCT) method. 

Area coverage and number density 
were obtained using the following meth- 
ods. We first isolated intergranular 
space using part of the procedure devel- 
oped by Crockett et al. (2009). While 
Crocketts procedure is optimal for some 
telescopes, the higher resolution of the 
NST required us to use a modified ver- 
sion method of Sanchez Almeida et al. 
(2004). Instead of playing series back 
and forth to single out the BPs we used 
NAVE method to track the plasma flow 
of individual BPs that could be tracked 
for at least 2 frames and were at least 
0.13". 

We used the wavelet technique to 
obtain the feature size power spec- 
trum. The wavelet technique was ap- 
plied following the constraints described 
in Andic (2010), however, we increased 
the confidence level to 99% to reduce 
the influence of noise induced irregularly 
across the FOV by the PCO.2000 cam- 
era, most prominently in the upper left 
half of FOV. 

Wavelet analysis is a very efficient 



filter;it separates different period oscilla- 
tory signals from processed signal (Tor- 
rence and Compo, 1998) forming a set 
of sub-ranges of oscillations across the 
main period range. For our spatial anal- 
ysis the main range is 0.076" - 11", and 
analysis formed 60 sub-ranges. Only pe- 
riodic signals that match set conditions 
are kept. Wavelet automatic analysis 
calculates confidence levels for every de- 
tected oscillatory signal. False detec- 
tions are reduced to the minimum when 
confidence levels are set to 99%. The fi- 
nite nature of the signal will cause false 
detection near the edge of signals. These 
false detections can be removed by ig- 
noring all oscillatory detections near the 
edge of a signal, an effective technique 
for shorter period oscillatory signal, but 
not for the longer periods. Prolonging 
a signal is the only viable option for 
the longer periods. We chose the size 
of sub-areas to 300 x 300 pixels as op- 
timal for our study, because we can no- 
tice the change in granulation shape and 
false detections are prominent only for 
structures of > 3" . 

We selected a segment of the FOV 
(Fig. [2]A) and analyzed 112 BPs over 
time. A segment was selected to avoid 
camera-induced noise, most prominent 
in the upper left portion of the FOV. 
Withing this segment all BPs who ap- 
peared during the time sequence and 
matched set conditions (size > 0.13" and 
time of the detection > 30 sec) were an- 
alyzed. Camera noise resulted in numer- 
ous artifacts of 0.1", which would un- 
dermined the reliability of our analysis. 
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Displacement of BPs 




arcsec 



Figure 2: An individual image of the selected region for analysis (panel A) 
and displacement of the analyzed BPs during their individual lifetime (panel 
B). Over granulation image the traces of BPs movement in time were plotted 
(blue solid line). Persistent conglomerations (i.e., knots) of BPs are visible 
at coordinates (8",9"), (25",6"), and (45",4"). 



3 Results 

3.1 Spatial distribution of 
BPs 

The spatial distribution of BPs across 
our FOV was analyzed following 
methodology presented in Sanchez 
Almeida et al. (2004, 2010) in a frame 



that had the highest difference between 
intensity of the bright and dark features. 
Because we counted more BPs than re- 
ported in previous research (Sanchez 
Almeida et al. 2010), we repeated anal- 
ysis on 7 more frames to increase re- 
liability of our results. We noted on 
average 743.62 BPs per frame, yielding 
2.8 BP per Mm^. 
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Figure 3: Size distribution of analyzed BPs. Y axis presents percentage of 
BPs in logarithmic scale. The largest percentage of BPs has 0.22" and that 
value is also close to the average value of diameter of all measured BPs. 



Percent coverage of a FOV depends 
on BP size (Fig. E]). If we assume 0.22" 
as the average size, BPs cover 5.38% 
of the FOV. This spatial distribution 
is the largest observed so far (Sanchez 
Almeida et al. 2010). Sanchez Almeida 



et al. 2004 stated that BPs are ubiqui- 
tous in the photosphere of the QS. We 
cannot confirm this, since we did found 
several patches of the 3.7" x 3.7" QS 
where we could not resolve any small 
features (Fig. [3|). 



3.2 BPs concentrations 

A majority of BPs were part of persis- 
tent knots that remained stationary dur- 
ing our observations. We define "knots" 
as locations where persistent BP's are 
observed for the duration of the time se- 
ries (Fig. [2p). 

The change in size distribution of 
granules could point to an infiuence of 
a constant presence of a BPs concentra- 



tions. We choose 4 distinct areas of the 
same size (300 x 300 pixels) for wavelet 
analysis; the QS (Fig. [U^) and area with 
ribbon as reference areas (Fig. [5]), and 
two areas with knots. 

We noted a structure stable during 
the time series (Fig. [5]), with the shape 
similar to ribbons previously noted in 
moth areas (Berger et al. 2004). Our 
observations were preformed during the 
solar minimum, and there were no active 
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Figure 4: Analyzed frame from dataset. We increased visibility of BPs by 
multiplying their intensity by a normalized Gaussian function with half-width 
0.14". 
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regions present. Also, during our obser- 
vations, SOHO MDI was not providing 
data. For several days prior to and af- 
ter our observational run, hoverer, in the 
same region MDI magnetograms showed 
only "salt and pepper" flux distribution 

QS area (Fig. [T]A.) was expanded to 
300 X 300 pixels to minimize an edge ef- 
fects induced by wavelet analysis. This 
expansion included sample of BPs in the 
QS. Nevertheless, even with additional 
BPs this reference area remained the 
quietest. 

Two remaining areas were 300 x 300 
pixels areas, each containing knots from 

The modal granular size was QS is 

I. 7" , with the smallest observed struc- 
ture as 0.2" (Fig. [6l solid black line). 
The sub-area containing a knot had 

II. 32% fewer structures in 0.8" - 2.7" 
than QS and the mode was 1.3" (Fig. [6l 
blue dashed line) . The sub-area contain- 
ing the second knot (Fig. El red dotted 
line) contained 15.48% fewer structures 
in 0.8" - 2.7" than QS and the mode was 
1.8". 

An area with a ribbon has 1.6" gran- 
ule as the mode with an abundance 
of structures < 0.8" (Fig. [6l orange 
dot-dashed curve). This area contained 
26.46% fewer structures in 0.8" - 2.7" 
than QS. The number of BPs per Mm^ 
in this area is 8.59. The size distribution 
curve for the area shows a noticeable in- 
crease in power for structures 3". This 
increase is in part caused by the dimen- 



typical for the QS regions. One observed 
ribbon and lack of magnetic information 
makes our dataset inappropriate for an- 
alyzing in detail this structure. We used 
this region as the control region, since it 
contained a plethora of BPs. 

coordinates: (8",9") and (25",6"), re- 
spectively. 

We analyzed the size of structures 
in each sub-area using wavelet analysis, 
yielding distribution of sizes (Fig. E]). 
Analysis was done separately for each 
frame in the time series and results are 
integrated over the time. 



sions and shape of the ribbon in addition 
to edge effects. 

Influence of conglomerations of BPs 
is most prominent in area with a rib- 
bon; which had 30.12% more structures 
in 0.12" - 0.56" than QS, while areas 
with knots had 3.54% more structures 
in 0.12" - 0.56" than QS. 

The dataset limits the possibility of 
connecting locations of knots to a net- 
work or intranetwork, because the posi- 
tion of network in the FOV cannot be 
determined precisely (Fig. S]). In ar- 
eas with knots, there were 3.24 BPs per 
Mm^ (Fig. m blue dashed line) and 3.29 
BPs per Mm^ (Fig. El red dotted line) 
making it 16% greater than average for 
the whole FOV. 
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Figure 5: A zoomed out right top corned of FOV. This structure is called 
a ribbon and it is considered to be conglomerate of the BP in areas with 
strong magnetic flux. With NST resolution we can resolve individual BPs 
and clusters of BPs that compose this ribbon. 



3.3 BP statistic 

A majority of analyzed BPs were in 
0.13" - 0.48" (Fig. ED. Larger BPs are 
rare; only 6% of all analyzed BPs were 
larger than 0.48". The mode was 0.22" 
comprising 19% BPs, while the average 
diameter of BP was 0.22". 

The lifetime of BP is influenced by 
abilities of instruments. Because BPs 
sizes are strongly affected by instrumen- 
tal resolution, measured lifetime is ac- 
tually the product of BP visibility to 



instruments rather than the period be- 
tween their physical creation and disap- 
pearance. BPs average lifetime is 4.28 
min, while median lifetime was 2 min, 
in a range of our imposed minimum of 
30 s to the longest lifetime of 29.75 min, 
while the modal lifetime is 0.5 min (Fig. 
[7|). The large percentage of BPs with 
lifetime of 30 s is likely a consequence of 
an imposed minimum. We also detected 
2% of BPs that were traceable trough 
their plasma flow for 30 min. All long- 
lived BPs were located in the knots. 
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Figure 6: Result of wavelet analysis of sub- images. Black solid line represents 
distribution of different sized structures in the granular subfield of QS. Blue 
dashed line represents same for sub-field that contains one of the analyzed 
knots, while dotted red line other knot. Orange dot-dashed line represents 
distribution of structures in subfield containing the ribbon. 



Intensity of each frame was normal- 
ized to the mean intensity of a flat field 
frame. We measured intensity of QS 
(Fig. [8l black line) and compared it with 
the intensity of BPs. (FiglHl blue line). 
Smaller areas of BPs caused larger er- 
rors during the intensity measurements. 



The QS intensity curve with its errors 
lie inside the error range for the inten- 
sity of the BPs. Thus, we can state that 
brightness of the BPs and QS are equal 
at the formation level for the continuum 
of the spectral line TiO 705.68 nm. 



Contrast averaged 0.1% for BPs, 
however, intensity error for the data-set 
was 0.07%, hence we can state that BP 



contrast in the continuum of the spec- 
tral line TiO 705.68 nm is equal to the 
QS contrast. 
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Figure 7: Distribution of visibility time of BPs within the dataset. Y axis 
presents the percentage of the analyzed BPs in logarithmic scale, x axis the 
time of the visibility. The visibility time of BPs do not match the true lifetime 
of the structures. 



4 Discussion 

4.1 Spatial distribution of 
BPs 

We used a broadband filter centered at 
TiO 705.68 nm spectral line, previously 
used only for an umbral investigations 
because it is very weak above granu- 
lation. Although use of TiO spectral 
line and broadband filter provide images 
consisting mainly of line continuum, our 
results can provide a reliable estimate of 
the distribution of small scale structures 
though not necessarily as a proxy for in- 
tense magnetic concentrations (Shelyag 
et al. 2004, Beck et al. 2007). No 
detailed analysis of the connection be- 
tween magnetic properties of small scale 
structures in TiO line has been done to 
our knowledge. 



A BP > 0.13" should appear in at 
least 2 sequential frames to be consid- 
ered for our analysis. NST has a clear 
aperture of 1.6m, achieving resolution 
close to 0.1" in the TiO 705.68 nm spec- 
tral line. However, the camera used 
and the speckle reconstruction method 
induced artefacts of a size approaching 
the diffraction limit, hence we ignored 
all structures smaller than 0.13". We 
observed 2.8 BPs per Mm^, ~ 3 times 
larger than the previous estimate for the 
G-band spectral line (Sanchez Almeida 
et al. 2010). Since we used very sim- 
ilar method as one used by Sanchez 
Almeida et al. (2004, 2010), we can 
speculate that differences in our results 
were caused by the different spectral line 
used. 
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Figure 8: Intensities of QS and BPs, normalized on a mean intensity of flat 
field frame. Black line with green error bars represent intensity of QS with 
corresponding errors. Blue line with gray error bars represent intensity of 
BPs. The intensity of BPs was obtained by measuring intensity changes of 
each BP and then averaging them. Intensity error is calculated using IDL 
procedure for the subsequent pixels. Error is equal to values where standard 
deviation of convolved area was smaller than standard deviation of same 
area convolved with a transposed kernel. These intensity profiles indicate 
that at TiO continuum formation level, intensities of BPs and granulation 
are similar. 



4.2 BP concentrations 

BPs showed a tendency to group at sev- 
eral locations in the FOV, which we call 
"knots". We analyzed 2 persistent BPs 
knots at the center of the solar disk (Fig. 
[2]). Knots seem to influence granula- 
tion in their immediate vicinity by caus- 
ing appearance of smaller granules (Fig. 
[6]). Power spectrum of structure sizes 
in 4 analyzed sub-areas showed a drop 
in power for structures in 0.8" - 2.7" in 
sub-area that with high concentrations 
of BPs. 

In subfields around knots the num- 



ber density of BPs was 16% larger than 
average for the whole FOV. In both ar- 
eas with knots, BPs cover 0.5% and 
0.6% more than average, indicating that 
mere spatial distribution of BPs cannot 
be the cause of the different size distri- 
bution (Fig. [6]) in those areas; rather 
their clustering around the knot loca- 
tion caused drop in power for structures 
in 0.8" - 2.7". On the other hand, the 
area with a ribbon has 8.59 BPs per 
Mm'^ and BPs cover 12.8% of sub-area 
(11.1" X 11.1"). In this case the shape of 
the size distribution curve (orange dot- 
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dashed line, Figl6]) is influenced by BP 
numbers too and not only by flux con- 
centrations. We can conclude that size 
of granules appears to be affected by the 
presence of a concentration of BPs in 
QS. 

Suppression of convection by mag- 
netic fields was known for decades 
(Parker, 1978). Morinaga et al. (2008) 
established that convection is sup- 
pressed by high concentrations of mag- 
netic flux tubes, not by strength of mag- 
netic fields. Arguments suggest that 
some kG flux tubes are not bright in the 
G-band (Sanchez Almeida et al. 2001, 
Vogler et al. 2005, Beck et al. 2007). 
This agrees with Ishikawa et al. (2007), 
who found areas of the high magnetic 
field without BPs. These arguments 
may indicate that higher resolutions are 
capable of resolving more kG flux tubes 
as BPs. 

If we assume that BPs in TiO con- 
tinuum are proxies for intense magnetic 
concentrations, we can speculate that 
our result is also based on flux tube con- 
centration; i.e. the noted drop in the 
granulation size is probably caused by 
the concentration of the flux tubes in the 
knot. 

This suppressing is localized to the 
vicinity of a BP knot. This effect is eas- 
ily noticeable and measured in active re- 
gions and network patches, since BPs 
are found over many arcseconds. It is 
more complicated to detect this effect in 
the QS since the cluster dimensions are 
usually comparable to the resolution of 
the instrument. With our increased res- 
olution we were able to measure those 



effects. 

4.3 BP statistic 

Crockett et al. (2010) stated that mag- 
netic BPs cannot be generated in large 
diameter magnetic flux tubes, and Beck 
et al. (2007) concluded that very few 
BPs are larger than 0.4". This is con- 
sistent with our observations that only 
6% BPs were larger than 0.5". With the 
resolution achieved in this dataset, we 
cannot see the real shape of an average 
BP, since the smallest resolved BPs are 
equal to our set limit. 

Viticchie et al. (2010) put forward 
the upper limit of ~ 0.1" as the dimen- 
sion over which the elementary G-band 
bright features are formed. With our ob- 
servations we could neither confirm nor 
dispute this. The average size of BPs 
in our dataset is two times larger, and 
the smallest BPs we observed was larger 
than limit put forward by Viticchie et al. 
(2010). 

The size of BPs in this dataset agrees 
with the result of Utz et al. (2009) and 
Beck et al. (2007). Utz et al. (2009) ob- 
tained different mean values for different 
sampling sizes (with sampling of 0.054 
arcsec/pixel authors measured 166 ± 31 
km). This agrees with our result of 
BPs diameter of 0.22" measured with 
sampling of 0.037 arcsec/pixel. Beck et 
al. (2007) observed BPs with an aver- 
age diameter comparable with our re- 
sult; differences in numbers were most 
likely caused by different sampling sizes. 

Rimmele (2004) observed 30%-50% 
larger apparent size of BPs in the G- 
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band than in the continuum. Although 
we analyzed BPs near continuum, we 
cannot confirm this result. We found 
that mean diameter of BPs is compara- 
ble with the size of BPs in the G-band 
reported by Rimmele (2004). 

Berger et al. (2004) observed larger, 
amorphous ribbons in an active region 
plage near disk center. We observed 
a similar structure in the QS. Unfortu- 
nately, at the day of our observations 
there was no magnetogram obtained by 
MDI, SOHO. The magnetograms from 
the previous and following days show 
the magnetic flux distribution typical for 
QS in the area where the ribbon was lo- 
cated. The difference between our struc- 
ture (Fig. [5|) and the one analyzed by 
Berger et al. (2004) is that we could 
see our ribbon resolved to the individ- 
ual BPs and clusters and ours was lo- 
cated in the QS. However, from analysis 
done with this dataset we cannot offer 
any explanation of this structure. Fu- 
ture analysis may benefit from inclusion 
of magnetic and Doppler-velocity infor- 
mation. 

The lifetime of BPs analysed is 
strongly influenced by the resolution of 
the instrument and the duration of the 
dataset. We had to impose lower tem- 
poral (30 s), and spatial (0.13") lim- 
its on the observed structures, limiting 
what we were able to see. The dura- 
tion between physical appearance and 
disappearance of BPs will be possible 
to establish only when spatial resolu- 
tion of instrument is high enough to re- 
solve an individual flux tube and follow 
it from formation to disappearance. In 



our dataset, the moment when BPs ap- 
pear or disappear does not correspond 
to BPs actual formation or disappear- 
ance of BPs, but is instead a measure of 
the visibility of BPs within the limits of 
our instruments. 

4.4 Irradiance 

Variation in total solar irradiance is 
caused by magnetic fields at the solar 
surface (Chapman et al. 1996, Lean et 
al. 1998, Filgge et al. 1998). Krivova 
et al. (2003) argued from their model 
that solar surface magnetism is respon- 
sible for solar irradiance changes. This 
argument can be tested by measuring 
contrast of BPs, pronounced in molec- 
ular lines along the blue side of the vis- 
ible spectrum and diminishing as one 
goes toward the red side of the spec- 
trum. At near-infrared wavelengths con- 
trast reaches negative values (Tritschler 
and Uitenbroek, 2006). According to 
Planck's law, for a deep red spectral 
line, a mean contrast of BPs at the solar 
disk center should be very close in value 
to the contrast of the QS. Tritschler 
and Uitenbroek (2006) presented con- 
trast measurements for several spectral 
lines that followed the prediction from 
Planck's law. 

The broadband TiO filter allowed us 
to access a t^qq = 1 level, where we can 
expect hotter granules than in G-band. 
Hence, contrast of BPs observed with 
this filter should be closer to the con- 
trast of the granules (Fig. [8|). 

The frequency range of the broad- 
band TiO filter falls into the deep red 
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part of the visible spectrTim. Tritschler 
and Uitenbroek (2006) analysed spec- 
tral lines that fall around the deep red 
part of the visible spectrum. Consider- 
ing the scale Tritschler and Uitenbroek 
(2006) established, our finding that the 
contrast of BPs is equal to the contrast 
of the QS, is in agreement with their re- 
sults. Hence, there is a possibility that 
even BPs in the quiet Sun make a con- 
tribution to the global solar irradiance. 

5 Summary 

Small concentrations of BPs seem to 
cause a reaction of convection in the 
quiet Sun, slightly suppressing convec- 
tion in their immediate vicinity. 

The NST at BBSO resolved a 
plethora of BP-like structures in the QS, 
around 2.8 BPs per Mm^ with TiO spec- 
tral line. Those structures tend to form 
persistent knots, which in turn forced 
the shape of granulation pattern to ad- 
just themselves around the knot loca- 
tion over the duration of the observa- 
tional run. Size distribution of granules 
in the immediate vicinity of knots points 
to a slightly suppressed convection due 
to the increased concentration of BPs. 

On average, analyzed BPs arc 0.22" 
in size and detectable for 4.28 min. 
Their averaged contrast is the same as 
the contrast of the QS. 

Thanks are due to the anony- 
mous referee whose comments helped 
improving this work. A A thanks 
J. T. Villepique for help with En- 
glish. We gratefully acknowledge sup- 



port of NSF (ATM-O745744 and ATM- 
0847126), NASA (NNX08BA22G) and 
AFOSR (FA9550-09-1-0655). 
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